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1.  Introduction 


The  problem  of  electromagnetic  (EM)  wave  scattering  for  waves  propagating  through  a 
plasma  sheath  surrounding  reentry  vehicles  and  vehicles  traveling  at  hypersonic 
velocities  at  high  altitudes  has  attracted  the  attention  of  many  researchers  [1-4].  In  this 
report  the  influence  of  low  frequency  wave  turbulence  excited  in  a  compressible 
supersonic  plasma  flow  with  velocity  shear  around  a  vehicle  is  analyzed.  A  system  of 
nonlinear  equations  consisting  of  the  momentum  equations  for  ions,  the  mass 
conservation  equations  for  ions  and  electrons,  and  an  equation  for  the  density  of  adiabatic 
electrons  is  presented.  These  equations  are  complemented  by  the  Poisson  equation  for  the 
electrostatic  potential  associated  with  the  excitation  of  ion  acoustic  perturbations  due  to 
the  presence  of  the  flow  shear.  Our  equations  appropriately  take  into  consideration  the 
effects  of  ion-neutral  collisions  and  ion  viscosity.  This  system  of  equations  is  solved 
numerically  and  wave  spectra  of  excited  low  frequency  perturbations  are  obtained.  We 
investigate  the  interaction  of  the  incident  EM  wave  and  low  frequency  turbulent 
pulsations  inside  a  plasma  sheath,  considering  them  relatively  small.  In  this  case  we  can 
separate  incident  and  scattered  waves.  At  the  same  time  characterizing  the  interaction  of 
the  incident  wave  with  the  turbulent  plasma  pulsations  requires  the  use  of  nonlinear 
equations  [5,  6].  The  nonlinear  current  generated  by  this  process  in  the  case  of  a  high 
frequency  incident  EM  wave  contains  only  the  electron  component.  On  analyzing  the 
dispersion  equation  for  electromagnetic,  electrostatic  Langmuir,  and  ion  acoustic  waves 
involved  in  the  scattering  process,  we  find  that  different  types  of  scattering  scenarios  are 
possible.  In  the  case  when  the  frequency  of  the  EM  signal  is  much  greater  than  the 
electron  plasma  frequency  ( a>0  »  cope ),  the  scattered  waves  will  also  be  electromagnetic 

with  frequency  co0  ±  coA ,  where  coA  is  the  frequency  of  the  ion  acoustic  wave  generated  by 
the  excited  turbulent  wave  spectrum.  In  the  case  when  the  frequency  of  the  incident  EM 
signal  is  close  to  the  electron  plasma  frequency  ( co0  >  cope ),  two  processes  can  take  place: 

1.  Scattering  of  the  incident  EM  signal  into  another  EM  wave  with  the  frequency 


2.  Transformation  of  the  incident  EM  signal  into  a  Langmuir  wave  with  frequencies 
in  the  range  co0±coA. 

We  find  that  turbulence  caused  by  such  ion  acoustic  wave  instabilities  can  strongly 
influence  the  properties  of  electromagnetic  waves,  potentially  impacting  the  quality  of 
communication  channels. 

The  following  is  the  layout  of  this  report.  In  Section  II,  we  present  the  multi-fluid 
equations  underlying  our  linear  stability  analysis  and  the  nonlinear  saturation  stage  of  the 
instability  associated  with  a  velocity  shear  in  the  plasma  flow.  In  Section  III  we  discuss 
numerical  solutions  of  the  nonlinear  equations  including  a  description  of  the  spectra  of 
generated  density  perturbations  associated  with  the  ion  acoustic  waves.  In  Section  IV  we 
analyze  the  scattering  of  the  incident  EM  wave  on  the  electron  density  perturbations  and 
calculate  the  amplitudes  of  the  scattered  waves  and  scattering  cross  sections  for  the  two 
cases  described  above.  We  conclude  the  report  with  a  summary  and  a  discussion  of  the 
main  results. 

2.  Instability  of  a  plasma  flow  with  a  velocity  shear 

In  the  present  section  we  investigate  the  excitation  of  vortices  and  ion  acoustic  type 
perturbations  in  a  compressible  supersonic  plasma  flow  with  velocity  shear.  A  system  of 
nonlinear  equations  consisting  of  the  momentum  equations  for  ions,  the  mass 
conservation  equations  for  ions  and  electrons,  and  an  equation  for  the  density  of  adiabatic 
electrons  is  presented  below.  These  equations  are  complemented  by  the  Poisson  equation 
for  the  electrostatic  potential  associated  with  the  excitation  of  ion  acoustic  perturbations 
due  to  the  presence  of  the  flow  shear.  Equations  presented  below  also  properly  take  into 
consideration  the  effects  of  ion-neutral  collisions  and  ion  viscosity.  We  are  interested  in 
the  excitation  of  turbulent  pulsations  in  a  compressible  supersonic  plasma  flow  with 
velocity  shear  in  the  two-dimensional  case  (see  Figure  1  below). 
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Fig.  1.  Plasma  flow  has  velocity  v0  (jc)  directed  along  the  y  axis  and  velocity  shear  is 
along  the  x  axis. 

The  basic  equations  which  describe  the  plasma  flow  are  as  follows: 

1.  Equation  of  motion  for  the  ions: 

=  ~ZenP  <P-Vpi+  m,  ni  VA  V,  -  K  ni  Via^i  - \)  (D 

2.  Equation  for  the  adiabatic  electrons: 

TVne  =eneV(p  (2) 

3.  Mass  conservation  equation  for  the  ions: 

dn 

— L  +  div(ni  V; )  =  0  (3) 

dt 

4.  Poisson  equation  for  the  electrostatic  potential: 

A (p  =  4 ne{ne  -  n( )  (4) 


d\ 
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In  equations  (l)-(4)  the  flow  velocity  can  be  presented  in  the  following  form: 

vi  =  %  0> +  S\i  (x,  y,t)  (5) 

In  equation  (5)  v0y(x)  is  the  initial  flow  velocity,  which  varies  along  the  direction 
perpendicular  to  the  flow  direction,  and  Svt(t,x,y)  are  the  perturbed  velocities  which 
appear  due  to  excitation  of  vortices  and  ion  acoustic  oscillations  inside  the  plasma  flow. 

To  solve  (l)-(4)  numerically  and  to  find  excited  density  and  electric  field  perturbations 
and  their  spectra  we  first  introduce  dimensionless  variables  and  rewrite  the  system  of 
equations  in  dimensionless  form.  For  the  ion  dimensional  mass  density  pi  and 

dimensional  number  density  AT  we  introduce  dimensionless  mass  pt  and  number 
densities  Ni  as 

pi(x,y)  =  p0mmpi(x,y)  N{  =n0max  N{(x,y)  p()mdX  =  M,n0max  (6) 

where  p0max  and  n0max  are  maximum  dimensional  mass  and  number  densities  and  M.  is 
the  dimensional  ion  mass. 

Now  we  introduce  dimensionless  time  according  to 

*  =  <V>  (7) 

where  the  dimensional  ion  plasma  frequency  is  given  by 


and  t  is  dimensional  time. 

We  normalize  the  space  variables  as 

x  =  Atx  y  =  A .y  ,  (9) 

where  A(  is  the  distance  between  two  grid  points. 

For  the  dimensionless  ion  velocity  in  the  flow  we  use 

U;  (x,  y)=cbpAM.  (x,  y).  (10) 

For  the  dimensionless  electrostatic  potential  of  the  excited  perturbations  we  introduce 
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(11) 


MjftrA2 

0(x,y)=  1  *  ;  0(x,y)  . 

Ze 

Finally,  for  the  dimensionless  mass  and  number  density  we  use 

-N  -N. 

/))  max  "n  nia\ 

Now  we  can  rewrite  the  main  system  of  Equations  14  in  the  dimensionless  and 

conservative  form  as  follows: 

Momentum  equations: 

+  —  (Nuixuix )  +  —  (A Tu  u  )  = 

3  X  X  i  ix  ix '  ^  x  i  ix  ty ' 

=  -N^.-!k^L-&LfiLN(U  ~U  )  +  -^N(—  +  —)U 

A.  ox  m.  cd  .  co  A.  ox  oy 


d(Ntuiy)  d 


-  H - (A ^u.u  ) h - (A I-u.u  )  = 

dt  dx  y  dy  y  y 


Pin  Yh 

dy  A2  dy  m,  cbDi 


d2  d2 


=  -*—■ -Yk^L-^Ni(uiy-uny)  +  ^J~yNi(-J  +  -J)uiy 


A,  dx  dy 


Mass  conservation  equation: 


dNt  d  d 

~yr~  +  ~d~  (Niuix  )  +  —(Niuiy)-Q 
dt  dx  dy 


Poisson  equation: 


AO  =  —  (Z  -  N.)  +  —  J0,\  ArO 
Z  Z(Te/M.) 


Equation  for  adiabatic  electrons: 


~  (  £0  |  ~  £0  ~  M  co 

Ne=N0e  exp \j- \-l[  =  N0eT  =  N0e^ 


M  co2  A2 

1  Pl  1  (J) 


In  Equations  13-17  we  also  used: 


N.=Zf), 

p  M, 
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where  Z  is  the  charge  number  of  ions. 

Using  the  above  dimensionless  system  of  equations  we  investigate  the  linear  and 
nonlinear  stages  of 

1.  Vortices  with  Vxv^O 

2.  Ion  acoustic  fluctuations. 

We  choose  the  following  Gaussian  initial  velocity  profile  (see  Figure  2): 

v0>,  (x)  =  exp  j-  —  ^  |  (19) 

where  xmid  =  \Lx  ,  xwid  =  |  Lx  ,  and  Lx  is  the  size  of  the  simulation  box  in  the  x 

direction.  Note  that  this  profile  has  inflection  points  (a  necessary  condition  for  the 
excitation  of  vortices)  and  it  can  also  be  responsible  for  the  excitation  of  ion  acoustic 
waves  inside  the  plasma  flow. 
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Fig.  2.  Gaussian  initial  velocity  profile  (red)  and  its  derivative  (blue)  for  the 
dimensionless  speed. 


Below  we  present  an  algorithm  for  the  numerical  solution  of  the  system  of  Equations  13- 
17  initialized  with  the  Gaussian  flow  velocity  profile  (see  Equation  19).  The 
dimensionless  time  evolution  equations  are  solved  using  temporal  and  spatial  finite 
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differences.  They  are  integrated  in  time  using  a  robust  two-step  predictor  corrector  with 
an  optional  Lax  step  in  the  predictor  [7,  8].  The  time  integration  scheme  is  illustrated  here 
using  the  density  equation,  where  the  density  is  denoted  by  N  and  the  superscripts  denote 
time  levels  as  n*At  (or  just  *  for  the  predictor  step  time  index)  and  the  subscripts  spatial 
locations  as  i  *  Ax  and  j*Ay,  where  A  s  indicate  the  unit  steps  in  time  and  space. 


Predictor: 


K,  =(v",)-(0.5*Ar)x{iL(ATl/ );+i  .  -(W,)’,  , 

-  (0.5  *  At)  x  {f  [( AV,)"  ,+1  -  (NVjI  M 

-  (0.5  *  A r)  X  D  X  { A£, -  2  *  AT  +  AT  w } 

-  (0.5  *  At)  x  D  x  { AT,,  -  2  *  N" .  +  AT,} 


Optional  Lax  Step  in  Predictor: 


{N"j)  =  + v",+1  +  n’h' 


No  Lax: 


Corrector: 


N"j'  =  V",  -  A(  x  j  i 


-(At)x{iL(wI);.ti-(jw<)*._i 
-(A()xDx{w”1J-2»Wj>iV"1J} 
-(A?)xZ>x{w,';M-2*W”J  +  Af“._1} 


(20) 


(21) 


Note  that  diffusion  acts  only  on  the  perturbed  part  of  the  density. 

One  unique  feature  of  our  numerical  algorithm  is  the  use  of  well-established  templates 
which  are  repeatedly  employed  to  perform  first  and  second  derivatives  in  space-centered 
finite  difference  form  in  both  the  x  and  y  directions.  These  are  the  only  derivatives  used, 


7 


given  the  conservative  formulation  of  the  convective  nonlinearities  in  all  three  dynamical 
equations.  Also,  the  Poisson  equation  is  solved  using  Fast  Fourier  Transforms  as 
incorporated  in  the  FFT2  set  of  routines  long  used  in  particle-in-cell  codes  developed 
over  many  decades  at  UCLA  [9].  Thus  the  Poisson  equation  reduces  to 

<22) 

with  nk  representing  the  perturbed  density  in  k-space  only  and  its  multiplier  constituting 

the  so-called  form  factor  where  an  appropriate  digital  filter  can  be  incorporated  if  needed. 
Furthermore,  a  dynamic  time  step  of  integration  has  been  implemented  in  our  algorithm 
whereby  the  instantaneous  time  step  is  adjusted  to  the  inverse  of  one-half  of  the 
maximum  x  or  y  velocity  anywhere  on  the  grid  thereby  insuring  that  the  Courant  stability 
and  accuracy  condition  [7,  8]  is  satisfied  at  all  time  in  the  nonlinear  calculations. 
Boundary  conditions  other  than  periodic  (and  even  periodic  for  ease  of  testing  of  the 
bounded  implementation)  have  been  incorporated  in  our  algorithm  using  the  image 
charge  method  [10].  This  means  that  the  appropriate  image  system  is  replicated  from  the 
physical  one  with  the  correct  parity.  As  an  illustration,  if  the  derivative  of  the  density  is  to 
go  to  zero  at  the  physical  boundary  or,  put  in  another  way,  if  the  density  is  continuous 
across  that  boundary,  then  the  density  in  the  image  system  is  such  that 

N(2Lx  -  x)  =  +N(x) .  (23) 

Periodic  boundary  conditions  imposed  over  the  extended  system  result  in 

N (Lx  +  x)  =  +N(x)  ,  (24) 

whereas  density  going  to  zero  at  the  wall  would  look  like: 

N(2Lx  -  x)  =  +N(x) .  (25) 

So  doing,  the  boundary  is  actually  between  the  last  grid  point  in  the  physical  system  and 
the  first  grid  point  in  the  image  system  taking  the  right-most  boundary  as  an  example. 
The  parity  for  the  other  terms  in  the  equations  consistently  follows  from  the  choice  made 
for  N.  Even  though  the  image  method  results  in  twice  the  expense,  it  is  much  to  be 
desired  for  symmetric  boundary  conditions  (i.e.  the  same  boundary  conditions  imposed  at 
both  boundaries)  because  of  its  simplicity  and  ease  of  implementation.  Results  below  are 
for  the  flow-driven  ion  acoustic  instability,  with  uniform  density  and  Gaussian  velocity 
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profile,  and  Te=T,=  1 .0.  The  parameters  are  such  that  VTi  —  copj  x  A; .  With  Te=Tj  it  also 

follows  that  the  unit  grid  spacing  A(.  =  rDe  where  rDe  is  the  Debye  radius  of  electrons. 

The  plots  below  (Figures  3-6)  are  for  maximum  flow  speed  equal  to  2.5  times  the  speed 
of  sound: 

W0max  =  2-5  X  Cs  .  (26) 

In  Figure  3  the  excitation  of  ion  acoustic  oscillations  by  a  flow  with  a  velocity  shear  in 
the  form  of  a  Gaussian  profile  is  presented.  The  black  curve  corresponds  to  the 
normalized  kinetic  energy  of  the  flow.  The  blue  and  red  curves  correspond  to 
dimensionless  density  and  electrostatic  potential  versus  dimensionless  time,  respectively. 


Fig.  3.  The  black  curve  corresponds  to  the  normalized  kinetic  energy  of  the  flow.  The 
blue  and  red  curves  correspond  to  dimensionless  density  and  electrostatic  potential  versus 
dimensionless  time,  respectively. 

It  is  clear  that  the  source  of  the  energy  for  the  instability  development  is  the  kinetic 
energy  of  the  plasma  flow.  Ion  density  and  electrostatic  potential  have  the  same  growth 
rate.  Saturation  of  instability  in  the  nonlinear  stage  is  connected  with  ion  dissipation. 

Growing  perturbations  of  the  electrostatic  potential  s>  at  different  moments  in  time  are 
presented  in  Figures  4  and  5.  According  to  Equation  17  the  electron  density  perturbation 
Ne  follows  the  perturbations  of  the  potential. 
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Fig.  4.  Potential  at  copit  =  900 
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Fig.  5.  Potential  at  copit  =  3000 
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The  formation  of  large  scale  structures  in  the  potential  can  be  observed  in  Figure  5.  Note 
that  these  density  structures  are  centered  in  the  regions  of  maximum  velocity  gradient, 
associated  with  the  Gaussian  velocity  profile.  In  Figure  6,  the  kx  and  ky  spectra  of  the 
density  fluctuations  are  presented.  The  decay  index  or  of  the  amplitude  of  the  density 
fluctuations  with  respect  to  normalized  kx  and  ky  is  found  to  be  within  the  interval 
-2. 15  <  a  <-1.98. 


Fig.  6.  Spectra  in  kx  and  ky  of  the  density  fluctuations 


The  results  of  the  spectrum  analysis  also  show  that  the  decay  index  of  the  ^-spectrum  is 
highly  dependent  on  the  temperature  ratio  %  and  on  the  maximum  of  the  flow  velocity 

^Omax  * 


3.  ELECTROMAGNETIC  WAVE  SCATTERING 

Below  we  present  an  analysis  of  the  scattering  of  EM  waves  on  density  perturbations 
associated  with  ion  acoustic  turbulence  excited  by  a  plasma  flow  with  velocity  shear.  The 
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amplitudes  and  spectra  of  the  scattered  waves,  as  well  as  the  scattering  cross  section  for 
the  EM  wave  propagating  through  the  plasma  sheath,  are  evaluated.  The  mechanism  of 
generation  of  scattered  waves  is  connected  with  appearance  of  induced  charges  and 
currents  under  the  action  of  the  incident  EM  wave.  We  will  investigate  the  interaction  of 
incident  wave  and  low  frequency  turbulent  pulsations  inside  a  plasma  sheath,  which  are 
assumed  to  be  relatively  small.  In  this  case  we  can  separate  the  incident  and  scattered 
waves.  On  the  other  hand,  the  process  of  the  interaction  of  the  incident  wave  with  the 
turbulent  plasma  pulsations  requires  the  use  of  nonlinear  equations.  The  nonlinear  current 
generated  by  this  process  in  the  case  of  a  high  frequency  incident  electromagnetic  wave 
contains  only  an  electron  component.  The  dispersion  properties  of  waves  involved  in  the 
scattering  process  are  presented  in  Figure  7,  where  the  wave  frequency  is  schematically 
plotted  as  a  function  of  the  wave  number  k .  It  follows  that  different  scenarios  are 
possible  from  the  analysis  of  the  dispersion  of  the  electromagnetic,  electrostatic 
Langmuir,  and  ion  acoustic  waves  involved  in  the  scattering  process. 
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Fig.  7.  Dispersion  curves  of  waves  inside  the  plasma  sheath: 

1  -  electromagnetic  wave  with  the  dispersion:  ^  _  ^J^+^c2 

2  -  Langmuir  wave  with  the  dispersion:  m  =  m  (\  +  -k2r2 )  »  where  r  =  Yu.  is  the 

/.  pe  2  De  De  CO 

pe 

electron  Debye  radius,  VTe  is  electron  thermal  velocity  and  cope  is  the  electron 
plasma  frequency. 

3  -  Ion  acoustic  wave  with  the  dispersion:  0)  _  %  where  co  .  is  the  ion  plasma 

frequency. 


In  the  case  when  the  frequency  of  the  incident  EM  wave  is  much  larger  than  the  electron 
plasma  frequency  (o0  »  a>pe ),  the  scattered  waves  will  also  be  electromagnetic  with 

frequencies  in  the  range  co0±a>A,  where  frequency  coA  belongs  to  the  wave  from  the  ion 
acoustic  wave  spectrum  excited  in  the  plasma  flow.  In  the  case  when  the  frequency  of  the 
incident  EM  wave  is  close  to  the  plasma  frequency  ( >  co  )  two  processes  can  take 

place: 
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1.  Scattering  of  the  incident  EM  wave  into  another  electromagnetic  wave  with 
frequency  oj()  ±  coA  . 

2.  Transformation  of  an  incident  EM  wave  into  a  Langmuir  wave  with  frequencies 
in  the  range  co0±coA- 


We  first  investigate  the  scattering  of  an  incident  EM  wave  with  frequency  a>0  »  cope .  We 
consider  a  plasma  sheath  with  the  following  parameters:  electron  density  n0e  =108  cm-3, 
electron  plasma  frequency  a>  =5.64 xlO8  rad  /sec,  ion  plasma  frequency 
copi  =  2. 1 1  x  106  rad  /  sec ,  and  the  frequency  of  the  incident  EM  wave  oj()  =  2/rx  1()9  rad  /  sec . 

The  density  of  the  neutral  air  molecules  inside  the  plasma  sheath  is  taken  to  be 
nN  =ioH  cm-3.  The  scattering  cross  section  of  electrons  on  the  neutral  air  molecules  for  the 

case  when  the  electron  energy  does  not  exceed  200  eV  can  be  taken  as 
cr  ~  1.5 xlO"15  cm2  [11].  For  a  temperature  of  electrons  Te  ~  0.01  eV  the  electron  -  neutral 

collision  frequency  can  be  estimated  as:  v  -nNVTae  =  6.29  x  102  sec-1.  In  this  report  we  focus 

attention  on  the  analysis  of  the  scattering  of  an  EM  wave  with  p  polarization.  This  wave 
can  be  presented  in  the  following  form: 

Eo  =  Eo,k0  exp{/(k0r  -  cokt) }  (27) 


where  for  the  components  of  the  electric  and  magnetic  fields  in  the  wave  we  have: 

Eo  ={£'o^£oy’°};  B0={0,Q,2U  .  (28) 


We  use  the  Maxwell’s  equations  to  describe  the  scattered  waves 

r,  1  dB+ 

c  dt 

471 .  1  SE+ 

rot  B.  =  —  i.  + - i 

*  c  ~  c  dt 

The  current  j±  contains  linear  and  nonlinear  parts: 


(29) 

(30) 


j± 


+J 


NL 


(31) 


In  order  to  find  the  linear  current  we  can  use  the  linearized  Vlasov  equation  for  the 
distribution  function  of  electrons: 


14 


dSf+  58  ff  e  1  dfM 

— — +  v — — - {E,  +-[vxB,])^!L  =  0  . 

dt  dr  m„  c  d\ 


In  Equation  32  Sf+  is  the  linear  part  of  the  distribution  function  of  plasma  electrons 
associated  with  the  electric  and  magnetic  field  of  the  scattered  wave.  fM  is  the  zeroth 


order  Maxwellian  distribution  function: 


/m  =  _ \3/2t73  eXP 


"  (2  kT1V^ 

where  u0  is  the  flow  speed. 


v?+(vy-«o,)2+v' 


The  linear  current  j+  is  the  first  moment  of  the  distribution  function: 


j±  =  ~ef  v8f+d} 


To  solve  Equation  32  we  will  choose  the  system  of  coordinates  moving  together  with  the 
plasma  flow: 

v,=vx;  V  =v y-u0  ;  VZ=VZ  (35) 


r  =  r  +  u0t-e 


where  v  and  r  are  in  the  new  coordinate  system  while  v  and  r  are  in  the  laboratory 
system.  In  the  new  system  of  coordinates  we  should  replace  the  frequency  of  the  incident 
wave  by  the  Doppler-shifted  frequency: 

O)0=a>0-  k 0yU0y  ’  (37) 

where  co  is  in  the  new  system  of  coordinates  and  d>  in  the  laboratory  system.  Now  it  is 
possible  to  write  equations  for  the  components  of  the  electric  field  in  the  scattered  waves 
in  the  following  form: 


On  introducing  6\(oj±,k±) ,  defined  as 


4-%-(^+^)=4^(^k+) . 

c  c  c 
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(41) 


we  represent  (38)  and  (39)  as: 


co+  .  .  .  _  .  A7tco+ 

-^e±(co±, k±)Ek±  =-i—^h±  . 

Now  we  calculate  the  nonlinear  current  j^L  which  is  responsible  for  the  scattered  waves 

generated  by  the  interaction  of  the  incident  EM  wave  with  low  frequency  ion  acoustic 
type  perturbations.  First  we  represent  the  distribution  function  F  of  electrons  in  the  form 

^  =  /*+<% +  «*+<%.  (42) 

where  fM  is  the  Maxwell  distribution  given  by  Equation  33,  Sf0  is  the  perturbed 
distribution  due  to  the  presence  of  the  EM  wave,  8f*  =  8f\±L  +  Sf\LN  are  the  linear  and 
nonlinear  parts  of  the  perturbed  distribution  function  due  to  the  presence  of  the  scattered 
signals,  and  SfA  corresponds  to  the  perturbation  of  the  distribution  function  due  to  the 

presence  of  ion  acoustic  turbulence.  The  equation  for  the  distribution  function  associated 
with  the  nonlinear  generation  of  scattered  waves  has  the  form 


d8f?  ,  „  88  j\ 


dt 


+  v  • 


1 


dSfA 


1, 


ddfo 


■  =  —  {E0  +  -[vxB0]}— ^  +  —  {£EA  +  -[vx5BJ} 
or  m  c  cy  m  c  o\ 


(43) 


In  Equation  43  E0  and  B0  are  the  wave  fields  of  the  incident  EM  wave,  SEA  is  the 
electric  field  in  electrostatic  ion  acoustic  perturbations,  and  SBA  =  0 .  We  use  Fourier 
transforms  to  solve  (43).  If  we  take  the  following  expression  for  the  incident  EM  wave: 


E0(r,  t)  =  Eko  sin(k0r  -  o\  t) , 
its  Fourier  components  are  given  as 

<Eo  >*,«=^:Ek/(ko  -k)^(<y0-<y)-^Ek  <J(k0  +k)^(ru0  +co) . 

2i  0  2 i  0 

The  Fourier  transform  for  the  ion  acoustic  waves  has  the  form 
8E7 


(44) 


(45) 


If  , 

= -  8E A(k,t)exp{i(a>At  -  k  r)d  r dt . 

(2  n) 


(46) 


The  Fourier  components  of  the  nonlinear  part  of  the  distribution  function  of  scattered 
waves  are  given  as 

ie  1  1  oSfA  ie  1 

- —^[(Eo+-VxBo)-^]  + - -  _  „  . 

m  co- k-v  c  oy  m  co- k-v  ov 


SfkN  = 

J  k,ft> 


88  f 


(47) 
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In  Equation  47  the  operator  Fk  a  stands  for  the  Fourier  transformation  applied  to  the 

expression  in  the  brackets.The  induced  plasma  current  due  to  the  interaction  of  the 
incident  EM  wave  and  ion  acoustic  waves  (responsible  for  excitation  of  scattered  waves) 
is  given  by 

=~eSy-  SfkNJ3y  •  (48) 


Using  Equation  47  for  Sj*  w  we  obtain 


S)^=Sil+Sj2+Sj 


33 


where 


*A=- 


ie2 


1 


me  co- k-vdv 

1  r  1 


f  E,  8f  .  ,d3k'da>' 

J  k0oo0  J  A, k-k  ,00-co 


(49) 


(50) 


Si2=-—\d\-v 
m  J 


and 


ie 

Sh= - v- 


co- k-v J  c 

1  8 


f-[vxB0,,  ,]  —  Sf  ,d3k'dco' ,  (51) 

J  n  L  o.k  O)  J  ^  J  A, k-k  ,co-a>'  ’  v  7 


me  ®-k-v  5v 
It  is  possible  to  show  that: 

I  c>  j2 1  k-v 


[fE,  .,,£>/  .  .d^'dco'] 

LJ  A,k  co  J  0,k-k  ,00-00'  J 


(52) 


\si  i  © 


«i 


(53) 


and 


Sh 


eA> 


me 


WL I  «  ml 


N 


(54) 


Taking  into  account  (53)  and  (54)  we  finally  have  for  the  current  which  is  responsible  for 
the  excitation  of  scattered  waves 

3/g2 
mea>' 

Now  with  the  help  of  Equation  41  we  have 

3  co2 


m  sn  J 


(55) 


e:  =- 


pe 


k  ,oo  t  2 

«0 


J"^0,k>  ^A,k-k>-ffl  ^  ^  ’ 


(56) 


where  E^  wis  the  Fourier  transform  of  the  scattered  electromagnetic  waves,  and 
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(57) 


,  <  k2c2 

6k ,a>  ~  1  2  2 

CO  CO 

The  spectrum  of  the  scattered  electromagnetic  field  (Equation  56)  in  fact  consists  of  two 
Ek_,«,_  and  Ek  w  ,  where 

k+  =k0  ±kA  (58) 

and  co±  =  co{)  ±  coA  (59) 

In  Equations  58  and  59,  k0  is  the  wave  vector  of  the  incident  EM  signal  and  k  ,  is  the 

wave  vector  from  the  turbulent  spectrum  of  the  ion  acoustic  waves.  Using  (56),  (58),  and 
(59)  we  obtain 

3  col 


T?  _  Pe 

_  t  2 


J^0,k \co'^n. 


A,k-k \ co-co' 


d  k'dco' 


S\,com  n0 

2(0 pe  Eq  A  A 

t  2  0  ■  L^U  k,,  <7,?k+k  J  . 

«0  2t 


(60) 


Instead  of  the  field  given  by  (60)  we  can  introduce  fields  E  k  m  and  E+k  a)  ,  where 
(k_,®_)  and  (k+,®+)  are  given  by  (58)  and  (59),  respectively.  Now,  taking  into  account 


t 

that  ~  ^k_  from  (60)  we  obtain 


*A 


2l  _,«®-  no 


A  similar  expression  can  be  obtained  for  E+k  .The  electric  field  in  real  space  is  given 
by 

E"  =  Jdk^d^E”  a  exp[/(k_r  -  cot) .  (61) 

From  Equation  61  it  is  clearly  seen  that  the  spectrum  of  scattered  waves  consists  of 
shifted  wave  vectors  and  frequencies,  proportional  to  the  wave  vectors  and  frequencies  in 
the  spectrum  of  ion  acoustic  waves  excited  by  a  flow  with  velocity  shear.  The  expression 
for  the  nonlinear  current  (Equation  55)  can  be  written  in  the  form 
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SiL  =  — |E  oX^nAt_k.„.d3k'dm' 

m  m  * 


meco 
'ico2 

Pe  EotH-i 


(62) 


8  7rn0o) 


k  0,®-<q) 


s'\ 


+\l0,(D+G)0 


]. 


Again,  from  Equations  58  and  59  we  get 


3  co 


ST». k  ..  =—^E0SnkA'WA 

67rnnco 


and 


(63) 

(64) 


$JV  =  j ^a^a^Jnm. .0-  exP[*(k  r  -  ®-0 

In  order  to  calculate  the  scattered  power  we  should  define  the  increase  in  the  energy  of 
scattered  waves  during  one  second  as 

SQ  =-Re\<SyN(r,t)-E  \r,t)>d3r  ,  (65) 

where  Re  means  the  real  part,  E  is  the  complex  conjugate  of  the  scattered 

electromagnetic  field,  and  <...>  denotes  time  average.  Further,  using  the  following 
relations, 

8)N(r,t)  =  ^^\dk_dco_  —  E0Sn\  k  exp[/(k_r-<yj)]  (66) 

8;rnA  J  ™  o  o  - 


3  cot 


and  E“(r,t)  =  — —  \dkdco 

8  ?rnn  J  _ 


co 


J _ 1_ 

co2  s’  . 

k_,<y_ 


■Eo Sn  k„-k_ 


exp[z'(k_r -co_t)\  (67) 


we  finally  obtain 


9tt32co4  .  , 

= - — £02  Re  f  dk_dco  _  — 

J  cl 


n 


t* 

S  .  CD 


r\Sn\  .  I2 

5  k0-k  _,<y0-6y  _ 


(68) 


with  the  density  spectrum  defined  as 


I  Sn 


\2=- 


1 


(ki+k^fcoy 

where  A  ~  0.1,  based  on  our  numerical  results. 

Knowledge  of  (68)  allows  us  to  calculate  the  scattering  cross  section  a  due  to  the 
presence  of  ion  acoustic  turbulence  as 

SQ~ 


<J  = 


(69) 


4  n 
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where  cr  is  defined  as  the  ratio  of  the  scattered  power  (defined  as  the  work  done  per  unit 
time  by  the  field  created  by  the  nonlinear  current)  to  the  energy  flux  of  the  incident  EM 
signal.  From  Equation  68  we  see  that  it  is  necessary  to  include  dissipative  terms  in  the 
expression  for  ,  otherwise  SQ  =  0  .  Taking  into  account  collisions  of  electrons  with 
neutrals  we  have 


co 2  k2c2  v  co 2 

pt  _  | _ pe_  _  L  ,  •  en  pe 

b  k,a>  L  2  2 

CO  CO 

6  l,k ,o)  1  2 


CO 


CO  CO 


V  CO 

_  •  en  pe 
6  2,k ,o)  1  3  ‘ 

CO 


Now  we  can  write 


Re(— )  =  -— 


2,k,o 


k,w  &  l,k,w)  +(£’2,k,©) 


(70) 


There  are  two  ways  to  calculate  the  integral  in  (68).  In  the  first  case,  which  we  can  call 
non-resonance,  scattered  waves  are  not  plasma  eigenmodes,  i.e., 

£^(co±,  k±)*0,  (71) 


and  this  results  in 


co±  ^  o\  . 

In  the  second  instance,  the  resonance  case,  we  have 

£-/(u4,k±)~0 


and  [5,  6] 


Re(- 


-)  =  -n 


S(co_-cok  ) 


k  ,a> 


ds[ 


dco_ 


co_=co{k_) 


(72) 

(73) 


(74) 


It  follows  from  the  dispersion  properties  of  the  EM  and  ion  acoustic  waves  that  (for  the 
range  of  frequencies  and  plasma  parameters  of  interest)  we  should  consider  the  non¬ 
resonance  case  (Equation  71).  Using  (70)  in  the  integral  for  the  scattering  cross  section 
given  in  (69)  and  performing  the  numerical  integration  one  obtains  the  following  results. 
When  the  width  of  the  plasma  sheath  is  varying  in  the  interval 

3  <  4  <  15  cm  (75) 

the  amplitude  of  the  scattered  waves  is  found  to  belong  to  the  interval 
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(76) 


10“5  <  —  <  10"1  , 

E0 

where  E0  is  the  amplitude  of  the  incident  EM  wave.  The  corresponding  scattering  cross 
sections  belong  to  the  interval: 

1(T17  <  cj±  <  10“8  cm2  (77) 

It  is  also  worth  mentioning  that  the  phase  change  in  the  scattered  waves  can  impact  the 
phase-locking-in  process  for  signals  from  GPS  satellites.  This  is  connected  with  the  fact 
that  the  absolute  values  of  the  wave  vectors  of  scattered  waves  k±  =k0  ±kA ,  due  to  the 

inequality  k0  <s:  k  , ,  can  be  responsible  for  a  significant  change  of  the  integrated  phase 

shift  y/  ~  |  k±dx . 

4.  TRANSFORMATION  TO  LANGMUIR  WAVES 

Below  we  present  an  analysis  of  the  transformation  of  the  incident  EM  signal  into  a 
Langmuir  wave  due  to  interaction  with  density  perturbations  associated  with  ion  acoustic 
turbulence  excited  by  a  plasma  flow  with  velocity  shear  around  a  hypersonic  vehicle.  The 
transformation  process  can  take  place  when  the  frequency  of  the  incident  EM  signal  is 
close  to  the  electron  plasma  frequency  inside  the  sheath.  This  corresponds  to  an  electron 
plasma  density  ne  ~  1010  cm-3 .  The  efficiency  of  the  transformation  of  an  EM  wave  (in 

our  case  it  is  the  GPS  signal)  propagating  through  a  plasma  sheath  into  an  electrostatic 
Langmuir  wave  on  ion  acoustic  density  perturbations  will  be  evaluated.  We  also  show 
that  inside  the  plasma  sheath  the  process  of  transformation  is  resonant  in  nature,  in 
contrast  to  the  EM  signal  scattering  into  an  electromagnetic  wave.  As  a  result  the  excited 
wave  amplitudes  as  well  as  the  scattering  cross  section  for  this  process  are  much  larger 
than  the  corresponding  quantities  in  the  case  of  scattering  into  an  electromagnetic  wave. 
The  mechanism  of  appearance  of  electrostatic  Langmuir  waves  is  connected  with 
transformation  of  an  incident  EM  wave  on  ion  acoustic  density  perturbations  by  induced 
charges  and  currents.  So  in  the  process  of  this  interaction  another  type  of  wave  -  an 

electrostatic  Langmuir  wave  -  is  generated  and  this  is  the  reason  it  is  called  a  wave 
transformation. 
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Dispersion  properties  of  waves  involved  in  the  scattering  process  are  presented  in  Figure 
1.  Below  we  concentrate  on  the  case  when  a>0  >  cope  (frequency  of  the  EM  wave  is  slightly 

above  the  plasma  frequency  inside  the  sheath)  and  transformation  of  the  EM  signal  into 
an  electrostatic  Langmuir  wave  on  ion  acoustic  density  perturbations  is  possible.  We  will 
consider  a  plasma  sheath  with  the  following  parameters:  electron  density 
n0e  ~  1010  cm-3,  so  that  the  electron  plasma  frequency  is  slightly  smaller  than  the 

frequency  of  the  incident  EM  wave  (GPS  signal  a>0  ~2ny\  09  rad /sec).  The  density  of 

neutral  air  molecules  inside  the  plasma  sheath  is  taken  to  be  nN  =10“  cm-3.  As  in  the 

previous  section,  we  analyze  the  scattering  of  an  electromagnetic  wave  with  p 
polarization  (see  (28)). 


The  scattered  Langmuir  wave  being  an  electrostatic  wave,  will  have  the  following 
polarization: 

Ej  ={Elx,Ely,0}\  ={0,0,0}.  (78) 

To  describe  the  scattered  wave  we  use  the  Poisson  equation: 

c//vE±  =  -Anerr ,  (79) 


where  E±=-VO±  and  AO1  =4 nerr  (80) 

The  density  rr  contains  linear  and  nonlinear  parts 

n±=nf+n*i.  (81) 

In  order  to  find  the  linear  density  perturbation  trL  we  use  the  linearized  Vlasov  equation 
for  the  distribution  function  of  the  electrons: 


dt  dr  me  dv 


(82) 


In  equation  (82)  Sf+  is  the  linear  part  of  the  distribution  function  of  plasma  electrons 


associated  with  the  electric  field  of  the  scattered  wave.  fM  denotes  the  zeroth-order 
Maxwellian  distribution  function  defined  as 
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(83) 


^=^rxp(- 


Vx+(Vy-M  0y)2+Vz 


where  u0y  is  the  flow  speed.  The  linear  density  perturbation  Snf  is  the  zero-th  order 
moment  of  the  distribution  function  given  as 

nj  =  J  Sf[d\  .  (84; 

To  solve  (79)  we,  as  in  the  previous  section,  choose  the  system  of  coordinates  moving 
together  with  the  plasma  flow,  and  hence  obtain 

Sft  =-—  k=,'  <V/„  (85: 

m„  ox  -k.-v 


<  = 


\SfLd’ 


en(A, 


k  \ 


(2 7r)my5Teme  J  m  -k±-v 


y2  +  y2  +  y2 

exp[ — \ - z-]d\  ,  (86) 

2vr, 


where  i\  represents  the  linear  part  of  the  density  fluctuations.  After  substituting  (86)  into 
(4)  we  finally  write  the  Poisson  equation  as 

sL  (k± ,  a>±  )£±2<£k  =  -Anenf  ,  (87) 


where 


with  w(z)  defined  as 


sL  (k± ,  co± )  =  1  +  — r— r-  [1  +  iyfirzwiz)] , 


w(z)=n^ptfi 

X  L  z-y 


yf2. k±yTe 


42yTe 


Finally,  the  Fourier  component  of  the  scattered  Langmuir  wave  field  can  be  written  as 


Ek  =-ik 


1  (k_-E0) 


_  _ _ PL  1  v  -  fi„  \r 

2  n0  G)2_  £L(k _,a>_)  k2_  Ka'°a  "  ' 


Hence  we  obtain  the  following  expression  for  the  scattering  cross  section: 
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36;r4A2  M 


°L  =■ 


dk^dk  do>t  C°S  6°  +  **’  S‘n  61(1)2  x 

Ir 


1 


(92) 


1^  t  f. 

\Z  KAx,KAy,(VA 


(O)0  -  COA )3  I  £\L  I2  +  I  ^2  I2 
Taking  into  account  that  scattered  Langmuir  waves  are  now  the  plasma  eigenmodes,  we 
make  use  of  expression  (74)  in  Equation  92.  Also  we  use  the  following  expression  for  the 
scattering  cross  section  for  the  numerical  integration: 


°Lres  = 


C  0)1 


s^sin^ 


kQ-k_,o)Q-(ok 


I2 


(93) 


with  the  density  spectrum  defined  as 
I  Snk  k  a  \2=^^~ 

KAx^KAy->a)A 


1 


ik^k2J  C0\  ’ 


(94) 


Ay  >  A 

where  according  to  numerical  results  A  ~  0.1.  Now,  when  the  width  of  the  plasma  sheath 
is  within  the  interval 


3  <  Lx  <  15  cm  (95) 

the  amplitudes  of  the  scattered  Langmuir  waves  are  found  to  be  much  larger  in 
comparison  with  the  amplitude  of  the  scattered  EM  waves.  The  ratio  of  the  scattered 
Langmuir  wave  amplitude  E'±  to  the  amplitude  E0  of  the  incident  EM  wave  was  found 
to  be  in  the  interval 

Fl 

10"2  <  -±-  <  1  .  (96) 

E 

Corresponding  scattering  cross  sections  are  in  the  interval 

10"5  <  cr*  <  10"3  cm2 .  (97) 
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5.  SUMMARY 


The  purpose  of  the  work  is  to  understand  the  influence  of  a  plasma  flow  with  velocity 
shear  on  the  properties  of  incident  EM  waves  used  for  communication  purposes.  This 
report  consists  of  two  parts.  In  the  first  part  we  have  investigated  the  excitation  of  low 
frequency  turbulence  inside  a  plasma  flow  with  velocity  shear.  We  have  demonstrated 
that  the  flow  with  a  Gaussian  velocity  profile  excites  low  frequency  turbulence  associated 
with  appearance  of  the  vortices  and  electrostatic  ion  acoustic  type  oscillations.  In  the 
second  part  of  the  report  we  analyzed  the  scattering  of  an  incident  EM  wave  on  excited 
turbulent  pulsations  associated  with  the  ion  acoustic  type  waves.  It  is  shown  that  in  the 
case  when  the  frequency  of  the  incident  EM  wave  is  much  greater  than  the  electron 
plasma  frequency,  the  scattered  waves  will  be  also  electromagnetic  with  frequencies  in 
the  range  o\  ±  coA ,  where  the  frequency  coA  belongs  to  the  wave  from  the  ion  acoustic 
wave  spectrum  excited  in  the  plasma  flow.  The  scattered  EM  waves  in  this  case  are  not 
plasma  eigenmodes  and  can  exist  only  as  forced  perturbations  in  the  vicinity  of  the  region 
where  they  were  generated.  It  is  also  worth  mentioning  that  a  similar  type  of  scattering 
(only  with  different  type  of  waves  involved)  was  observed  during  the  active  experiments 
in  the  ionospheric  plasma,  analyzed  in  [12,  13], 

In  the  case  when  the  frequency  of  the  incident  EM  wave  is  close  to  the  plasma  frequency 
( co0  >  cope )  two  processes  can  take  place: 

1.  Scattering  of  the  incident  EM  signal  into  an  electromagnetic  wave  with 
frequencies  a>0±coA. 

2.  Transformation  of  an  incident  EM  wave  into  Langmuir  waves  with  frequencies  in 
the  range  co0±coA- 

The  scattered  EM  waves  in  this  case  are  not  plasma  eigenmodes,  but  the  Langmuir  waves 
are.  As  a  result  the  scattering  cross  section  for  the  processes  involving  Langmuir  waves 
inside  a  plasma  sheath  is  much  higher. 
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